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The results of mathematical modeling of the thermal operation of a 280 tons/day glassmaking furnace with a 
horseshoe flame are presented. It is shown that specific extraction of molten glass 2.5 tons/(m 2 ■ day) and spe¬ 
cific heat flow 4563.8 kj/kg with gas heating of the furnace can be attained without additional electric heating 
of the melting tank. 


Key words: glassmaking furnace, mathematical modeling, heat-and-mass transfer, energy efficiency. 


The modem methods of designing glassmaking furnaces 
include mathematical modeling of the thermal work of a fur¬ 
nace. The difficulty in developing new furnace designs stems 
from the need to perform a detailed analysis of external and 
internal heat-and-mass transfer processes, which by their 
physical nature are not only coupled but also mutually de¬ 
pendent to a considerable degree [1]. 

Neither the heat balance of a furnace nor calculations of 
its individual structural parts reveal the coupled character of 
the thermal processes involved in glassmaking. Moreover, 
local methods of performing calculations do not permit de¬ 
termining the effect of the boundary conditions of forced 
convection of the melt on internal heat transfer and the struc¬ 
ture of the convection flows in the melting tank. The tempe¬ 
rature of the molten glass at the entrance into the neck of a 
glass furnace, which is the main boundary condition for cal¬ 
culating the extraction channel of furnace, remains undeter¬ 
mined. Of definite practical interest is the temperature distri¬ 
bution on the inner surfaces of the refractory masonry. These 
and other problems can be solved by mathematical modeling 
of the thennal operation of furnaces using coupled numerical 
models of the internal and external heat-and-mass transfer 
[2, 3], 

We shall examine the modeling results obtained for a 
glassmaking furnace with a horseshoe flame for melting con¬ 
tainer glass. The furnace design implements a number of 


1 First President of Russia B. N. El’tsin Urals Federal University, 
Ekaterinburg, Russia. 

2 E-mail: vdzuzer@yandex.ru. 


novel technical solutions. In combination they secure not 
only high glassmaking efficiency but also the possibility of 
making colored as well as colorless glass. In addition, even 
though the specific production 2.5 tons/(m 2 • day) of molten 
glass is high, the furnace design does not provide for addi¬ 
tional electric heating of the melting tank. 

The nominal energy efficiency of glassmaking 
(1090 kcal/kg) is attained by optimizing the organization of 
the flame and the high additional heating temperature of the 
air and minimizing the heat losses into the surrounding me¬ 
dium. The construction of the burner ensures that the condi¬ 
tions required for directed forward radiant heat transfer and 
fuel burnup in the range 70 - 80% of the length of the work¬ 
ing space of the furnace are created. 

Efficient heat-insulating materials are used for the crown 
and walls of a once-through regenerator. In consequence, the 
heat losses through the regenerator’s masonry do not exceed 
5% of the heat content of the products of combustion. The in- 



Fig. 1. Diagram of melting tank. 
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Fig. 3. Temperature field of the surface of the molten glass. 


Fig. 2. Temperature field of the surface of the roof. 


ternal volume of the regenerator is filled with TL 14/15 
pot-type elements. The specific surface area of heating of the 
charge is 45 m 2 /1 m 2 of the melting tank. As a result, the 
temperature of the exhaust gases at the exit from the checker 
is about 300°C, while the nominal amount of additional heat¬ 
ing of the air for combustion is greater than 1300°C. 

Technical Characteristics of a Furnace 


Capacity, tons/day (kg/sec). 280(3.24) 

Specific molten-glass output, tons/(m 2 • day).2.5 

Fuel.Natural gas 

Lowest working heat-generation capacity 

of fuel, kJ/m 3 . 33,496 

Air heating temperature in regenerator, °C.1300 

Gas flow to melting, m 3 /h.1590 

Specific heat flow to glassmaking, kj/kg. 4563.8 

Maximum furnace roof temperature, °C. 1580 ±10 


The results obtained in [5] were used to design the fur¬ 
nace weirs in the melting part of the furnace. The use of ener¬ 
gy-efficient masonry made it possible to reduce the heat 
losses into the ambient environment to 7% of the consump¬ 
tion part of the furnace heat balance. 

Special attention was devoted to the calculation and de¬ 
sign of the melting tank. The configuration solutions for the 
tank profile are based on previously established and mathe¬ 
matically formalized laws of hydrodynamics and internal 
heat transfer [6]. 

The scheme of the melting tank with melting and fining 
zones of variable depth separated by an overflow weir 
(Fig. 1) is well known in the glass industry and widely used 
in furnace design. At the same time its efficiency manifests 
only with a certain combination of geometric dimensions of 
the melting and fining zones, of which the most important 
are the height of the dam and neck as well as the length and 
depth of the pit after the dam. Of fundamental importance for 
forming a flow with a rational two-cycle structure in the tank 
is the degree to which these parameters match the boundary 
conditions for free (natural) and forced convection [1], 


t, °c 



Fig. 4. Variation of the surface temperature of molten glass aver¬ 
aged over the length of the furnace. 

A preliminary calculation of the melting tank is done 
using free-convection boundary conditions prescribed by the 
temperature distribution on the surface of the molten glass 
[7]. At this stage the depth of the melting and fining zones as 
well as the structure of the masonry of the side walls and bot¬ 
tom of the tank are calculated. 

In accordance with the concept of building a coupled 
mathematical model of the thermal operation of a glass¬ 
making furnace [1], the temperature field on the surface of 
the molten glass is determined by the external heat transfer in 
the working space. In the article the results of modeling the 
external heat transfer are presented by the temperature fields 
of the roof (Fig. 2) and molten glass (Fig. 3). 

Analysis of the data in Fig. 2 shows that the maximum 
heating temperature of the roof is 1582.2°C. This corre¬ 
sponds to the conditions of long-time operation of the heat- 
insulating silica roof [1]. 

The average temperature of the roof is 1485.8°C. 

The surface heating of the molten glass (see Fig. 3) is 
characterized by a relatively uniform temperature distribu¬ 
tion over the width of the tank in the melting zone. Even 
greater uniformity of the heating is observed in the fining 
zone (after the dam). The distribution along the furnace 
length of the tank-width average of the surface temperature 
of the molten glass (Fig. 4) corresponds to the conditions for 
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TABLE 1. Parameters of Heat Transmission through the Melting 
Tank of the Furnace (Dark-Green Glass) 



Melting zone 

Fining zone 

Parameter 

Free 

convection 

Forced 

convection 

Free 

convection 

Forced 

convection 

h, m 

1.2 

1.2 

2.0 

2.0 

: o/-t 

L mg » 

1373.7 

1373.7 

1391.1 

1391.1 

cr 1 

o 

o 

1305.2 

1317.2 

1279.4 

1338.9 

4xt, °C 

95.2 

95.7 

94.0 

96.7 

R, (m 2 • K)/W 

2.00 

1.89 

2.09 

1.88 

K, W/(m 2 • K) 

0.4784 

0.5058 

0.4596 

0.5085 

q, W/m 2 

638.0 

646.0 

621.0 

660.5 


Notations: h ) thickness of glass layer; f mg ) average temperature of 
the molten-glass surface; f b ) average temperature of the melt at tank 
bottom; f ext ) average temperature of the exterior surface of tank bot¬ 
tom; R) average coefficient of thermal resistance; K) average coeffi¬ 
cient of heat transfer through the masonry; q) average heat flux into 
the ambient enviromnent. 


the formation of two-cycle structure of the convection flows 
[1]. The temperature maximum lies at a distance from the en¬ 
try wall equal to 0.75 of the length of the melting tank. To fix 
the position of the hotspot along the tank a row of bubbling 
nozzles is placed precisely at this location. We note that the 
average (1380.4°C) and maximum local (1502.5°C) surface 
temperatures of the molten glass correspond to the condi¬ 
tions of high-temperature melting of glass. 

Analysis of the surface temperature field of the molten 
glass makes it possible to determine the boundary tempera¬ 
tures in order to calculate the melting tank. The average sur¬ 
face temperatures of the tank in the melting (before the dam) 
and fining zones are 1373.7 and 1391.1°C, respectively. Un¬ 
der these conditions and in the absence of glass extraction 
the melt temperatures at the tank bottom are 1305.2 and 
1279.4°C (Table 1), respectively. This temperature level at¬ 
tests to the high efficiency of the heat insulation of the bot¬ 
tom, which is confirmed by the very low heat losses into the 
ambient environment. For the melting and fining zone 
q = 638.0 and 621.0 W/m 2 , respectively. 

The real picture of the convective motion of melt in the 
tank can be obtained by taking account of the boundary con¬ 
ditions for forced convection. The streamlines, whose tan¬ 
gents indicate the direction of the velocity vector at the 
points of tangency, give a graphic picture of the instanta¬ 
neous flow pattern. For steady flow the streamlines are also 
the trajectories of the motion of the medium. 

The melt flow in the melting tank is represented by a 
two-dimensional field of relative flow lines (Fig. 5). The 
streamlines are normalized to the mass flow of the glass 
through the neck, which equals 3.24kg/sec for the proto¬ 
type. The magnitude of the relative streamlines characterizes 
the number of times the molten glass circulates in the corre¬ 



Fig. 5. Flow pattern of the molten glass in the melting tank (the 
stream functions are normalized to the flow rate of the melt through 
the neck). 


sponding rotational contours. It should be noted that there is 
a large difference of the intensity of the motion of the molten 
glass between the melting and fining zones. For this reason, 
the scale of the image for these zones of the melt tank is 
changed in order to clearly visualize the calculations. 

On the left-hand side of the tank (before the dam) there 
are two centers of circulation in the flow of the molten glass; 
the circulation multiplicity factor is 49.8 in the left-hand and 
16.54 in the right-hand center (it is not manifested in Fig. 5). 
On the whole the melt in the melting zone rotates counter¬ 
clockwise: from left-hand wall to the dam along the bottom 
of the tank and to the left-hand wall from the dam in the 
subsurface layer. 

Because the melt flow breaks off on the top edge of the 
weir a small contour of very intense circulation forms on 
left-hand surface of the dam with circulation multiplicity 
-82.2 at the center (the minus sign indicates clockwise rota¬ 
tion). Finally, after the dam there is a single center of very in¬ 
tense circulation with multiplicity 1633.9. 

The high circulation multiplicity of the molten glass after 
the dam promotes averaging of the melt temperature in the 
fining zone. In consequence, the coefficient of thermal uni¬ 
formity of the molten glass at the exit from the melting tank 
reaches 99.3%. The temperature of the glass at the entry into 
the dam is 1343.9°C. 

As a result of the intensification of convective heat trans¬ 
fer from the high-temperature fining zone into the batch 
loading zone the temperature of the melt layer at the bottom 
on this section of the tank (Fig. 6) is higher than at the sur¬ 
face of the tank (see Fig. 4). On the whole the average values 
of the bottom temperature of the molten glass in the melting 
and fining zones increase to 1317.2 and 1338.9°C, respec¬ 
tively (see Table 1). It should be noted that the temperature 
level indicated is reached, as a rule, with additional electric 
heating of the melting tank with power 1000 kW. Rational¬ 
izing the design of the melting tank and optimizing external 
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Fig. 6. Variation of the melt temperature at the bottom of the melt 
tank. 


heat transfer make it possible to eliminate losses due to the 
installation and operation of an auxiliary electric heater. 

In closing we can state that the mathematical modeling 
supports the validity of the design adopted for the furnace. It 
attests to the possibility of high-capacity energy-efficient 
glassmaking without an auxiliary electric heater. 
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